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Abstract: We study the evolution of curvature perturbations and the cosmic mi- 
crowave background (CMB) power spectrum in the presence of an hypothesized extra 
anisotropic stress which might arise, for example, from the dark radiation term in 
brane-world cosmology. We evolve the scalar modes of such perturbations before and 
after neutrino decoupling and analyze their effects on the CMB spectrum. A novel 
result of this work is that the cancellation of the neutrino and extra anisotropic stress 
could lead to a spectrum of residual curvature perturbations which is similar to the 
observed CMB power spectrum. This implies a possible additional consideration in 
the determination of cosmological parameters from the CMB analysis. 
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1. Introduction 

Anisotropic stress is the traceless component in the energy-momentum tensor and it 
accounts for anisotropic momentum flow in the universe. Although some previous 
works have considered neutrino anisotropic stress and the anisotropic stress of the 
primordial magnetic field Jj], |], || |], ||, ||, they mainly considered the post neutrino 
decoupling era. In this work we deduce the evolution of the neutrino anisotropic 
stress and curvature perturbations both before and after neutrino decoupling which 
can be important. 

In the standard cosmological model, the anisotropic stress is absent before the 
epoch of neutrino decoupling because rapid interactions among elementary particles 
dissipate it. After the neutrinos have decoupled from the cosmic expansion, how- 
ever, the neutrinos become relativistic free particles and neutrino anisotropic stress 
can grow gradually. The neutrino anisotropic stress plays an important role in the 
formation of large scale structure and the CMB fluctuation spectrum. 

In this paper, however, we consider other possible sources for anisotropic stress 
besides neutrinos or a magnetic field which can be present before neutrino decou- 
pling. In particular, such terms arise naturally in cosmological theories in higher 
dimensions. We follow the evolution of the neutrino anisotropic stress and curva- 
ture perturbations both before and after neutrino decoupling in the presence of such 
anisotropic stress terms and show that curvature perturbations can stay constant on 
superhorizon scales as does the standard adiabatic mode. However, unlike the stan- 
dard case, significant evolution of curvature perturbations occurs before neutrino 
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decoupling. We show that, for the right conditions, such curvature perturbations 
might even reproduce the observed CMB power spectrum implying a possible ad- 
ditional consideration in the extraction of cosmological parameters from the CMB 
analysis. 

One possible source of extra anisotropic stress is the so-called "dark radiation" 
term in brane-world cosmology. This term can affect strongly the CMB anisotropies 
[0, [8], |^, [TJJ, [IT], |T2j . In such brane-world cosmology, there are two correction terms in 
the 4D Einstein equation. One is the extrinsic curvature which introduces a quadratic 
term in the energy momentum tensor, but is negligible in the low energy limit. The 
other is from the projected Weyl tensor. The application of the five-dimensional 
conservation condition to the Weyl tensor requires that the energy density term 
varies with scale factor as a -4 to an observer on the brane Hence, it is called 
the dark radiation term and remains significant throughout the radiation-dominated 
epoch. 

There is, however, no intrinsic brane equation to determine the anisotropic stress 
||. Although there have been attempts (e.g. |TB| ) to constuct simple models for the 
Weyl anisotropic stress, the solution strongly depends upon unknown physics in the 
bulk dimension. Hence, in what follows we take a phenomenological approach and 
simply adopt the plausible assumption that in the same way that inflation generates 
scale-free fluctuations characterized by a spectral index P(k) ~ k n (with n ~ 1), 
we can expect that the Weyl anisotropic stress could be characterized by a spectral 
index to be determined (constrained) by a fit to the observed CMB power spectrum. 
Moreover, this dark-radiation term is assumed to scale as oc a -4 similar to the energy 
density and pressure as was adopted in Ref. [[|. 

Another more familiar example of extra anisotropic stress is that of a primordial 
magnetic field (PMF). The amplitude of the energy density B 2 /8ir and magnetic 
anisotropic stress p 7 7r# of the PMF again both scale as radiation density oc a -4 . 
Moreover, such a PMF should be characterized by an amplitude and spectral index, 
However, the contribution of a PMF to the CMB power spectrum is constrained to 
be small because it has a measurable effect on the observed BB mode of the CMB 



on small angular scales QT4J, [TJ|. As shown in Ref. |JLqj , the tensor mode of the 
primordial magnetic field has two kinds of fluctuations. One is a compensated mode, 
which arises from the compensation of anisotropic stress, and the other is a passive 
mode which was generated by the anisotropic stress of the magnetic field before the 
epoch of neutrino decoupling. The passive tensor mode makes a large contribution 
to the CMB. Although Ref. |16[ also studied the vector mode, there was no analysis 
of the scalar mode in that work. The spectral index of a PMF is also constrained 



from the effects of the associated gravity waves on primordial nucleosynthesis [17 



Hence, a PMF cannot source the scale-free extra anisotropic stress of interest to this 
work. We only mention it as an example of a power-law anisotropic stress which also 
scales as oc a~ A . 
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Here we show that a scale- free anisotropic stress, if present, could imply an 
additional consideration in the determination of cosmological parameters from the 
CMB. As an extreme example, in this paper we show that an extra anisotropic 
stress could even lead to a CMB spectrum that agrees with observation. Such a 
source is a bit contrived as it must be scale invariant as naturally occurs in the 
inflationary scenario. Thus some kind of mechanism would have to occur to produce 
the desired scale-invariant anisotropic stress. Our main point is that the possibility 
exists that such an anisotropic stress term might be present in the early universe and 
that it would lead to curvature perturbations which could reproduce the observed 
CMB power spectrum. This possibility has not previously been pointed out to our 
knowledge. 

A second point which we make here is that this extra source of anisotropic stress 
may produce non-Gaussian fluctuations depending upon the nature of the source of 
anisotropic stress. The WMAP-5yr analysis, has indicated that there is at least a 



possibility for non-Gaussianity although the uncertainty is large The Planck 

mission should constrain non-Gaussianity with high precision. If non-Gaussianity 
were actually observed, it could suggest (e.g. |L9]) the need for a new cosmological 
paradigm which allows for non-Gaussianity such as that described here. At the very 
least, observational limits on non-Gaussianity could place limits on the hypothesis 
proposed here, depending upon the source of the extra anisotropic stress. 

2. Evolution Equations 

To eliminate gauge freedom we study the evolution of the anisotropic stress in covari- 
ant coordinates within the frame of the cold dark matter (CDM frame). In this frame, 
the variables can be defined on a supersurface orthogonal to the CDM 4- velocity u^. 
Then, using the notation of Ref . || , the anisotropic expansion rate (shear) a and the 
inhomogeneity of the expansion Z can be determined from the covariant derivative 
of u^. We can neglect vorticity in the scalar mode. We also introduce the Weyl 
tensor, i.e. traceless part of the Riemann curvature tensor, which vanishes in the 
FRW background spacetime. The magnetic part of the Weyl tensor is negligible in 
the scalar mode. Therefore, we can define the electric part of the Weyl tensor as 
$. By linearizing the Bianchi and Ricci identities, one can obtain the following five 
equations for $, o and Z [g, [ZD| [HJ: 

$' + %® + —Kpa 2 (ja + i? 7 g 7 + R u q u + R b qb + -R 7 g ex ) 

npa 2 (37 — l)(i? 7 7r 7 + R u tc u + -R 7 7r cx ) 



2' 



2k 

-7^ K pa 2 (R^ + R u -k' u + Rj7r' cx ) = , (2.1) 
a' + Ha + k® + ^-Kpa 2 (R 1 Ti 1 + R u tc u + R^) = , (2.2) 
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Z' + HZ + -^-/tpa 2 (i? 7 A 7 + R V A U + R b A b + R C A C + i? 7 A cx 

+3{R y SP J + R v bP v + R b 5P b + RjSP ex )) = , (2.3) 
2 1 

-(Z - a) + — /«pa 2 (i? 7 g 7 + R u q u + R b q b + -R 7 g ex ) = , (2.4) 

2$ - —rKpa (.R 7 A 7 -|- R v Ay -\- R b A b -\- R C A C -\- R^A CX -\- /? 7 7r 7 -|- R u tc u -\- _R 7 7r ex ) 

— — npa 2 (R 1 q 1 + R u q u + R b q b + -R 7 g ex ) = . (2.5) 

Here, the prime denotes a derivative with respect to conformal time r and H = a'/ a. 
We have defined an equation of state parameter 7 as p = (7 — l)p, where p is the 
total pressure and p is total energy density, and the gravitational constant is written 
k = 8nG. The total energy density p is written p = p 7 + p u + p b + p c , and we 
have introduced the energy fractions R, e.g. i? 7 = p 7 /p, R v = Pu/p, Rb = Pb/p and 
Rc = Pc/ P- The subscript denotes each component, i.e. the photon (7), neutrino 
(u), baryon (b), CDM (c) and the extra source (ex) contributions. For normal fluid 
variables, the density fluctuations A, the heat flux q, and the pressure fluctuations 
5P are normalized using their background fluid energy density, e.g. p 7 A 7 = 5p 7 , 
p u A u = 5p u , etc. For the fluctuations related to the extra source A ex , <?ex and 7r cx , 
however, they enter the energy-momentum tensor directly without an associated 
background fluid energy density. Since we assume that they scale as a -4 , we can 
conveniently normalize to the photon energy density. Hence, in analogy with the 
other fluctuations, in the Fourier space, we define A ex , q ex , and 7r cx by 

p 7 A cx = -T CX ° (A:) , (2.6) 
p 7 g ex = k'Tjm , (2.7) 

P^ex = - -VX^/C) , (2.8) 

where ki is the wave vector with a unit length and T ex ^ is the energy-momentum 
tensor for the extra stress. Note that, because we write the extra ansiotropic stress 
as p 7 vr cx , p 7 carries the required a -4 scaling. Therefore, 7r ex becomes a fixed normal- 
ization constant. As noted above, the spectra of A ex , q ex and 7r ex are not given by 
the Einstein equation and require a specific model of the extra dimensions. 

Well into the radiation dominated epoch, the background energy densities for 
the baryons and CDM are negligible. Here, in this epoch, we take p ~ p 7 + p u , 
i? 7 + R u ~ 1 and the R 1 and R v are constant. Since we focus on the early universe, 
we need to study the evolution of perturbations on super horizon scales, i.e. kr <^ 1. 

The curvature perturbation r\ is defined as 

77 = -(2<5> + a'/k) . (2.9) 
7] also satisfies the relation 7/ = — \{Z — a) M. 
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Using the fact of npa 2 = 3H 2 and k(37 — 2) pa 2 = —6H', we can rewrite Eqs. 
ID- COD as follows: 



& + n$ + ^{U 2 - H')a + ^H 2 (i? 7 g 7 + R v q v + R b q b + R^q ex ) 



: 



(H 2 - 2U')(R 7 %~ ( + R„tt v + i? 7 7T cx! 



3T-L 2 

-^-(i? 7 7T 7 + R v -k' v + RjTc' cx ) 



3V, 2 

a' + Ho + &;$ + -^-(i^TT 



2A- 



7"7 



R V 1X V + -R 7 7T 

t 



■7 ' 1 ex ) 



m 2 



Z' + HZ + —(RyA^ + R U A U + R b A b + R C A C + P 7 A C 

+3( J R 7 5P 7 + R V 6P V + R b 5P b + R y 5P ex )) = , 
2 31-i 2 

-(Z - a) + -p-(i? 7 g 7 + + i2 6 g 6 + i? 7 g ex ) = , 



(2.10) 
(2.11) 

(2.12) 
(2.13) 



2$ 



3H 2 
k 2 



(i? 7 A 7 + R V A V + i?feA fc + R C A C + i? 7 A ex + _R 7 7r 7 + + -R 7 7r cx; 



97^ 

' A; 3 



(i? 7 g 7 + R u q u + -Rtgb + -R 7 g cx ) = . 



(2.14) 



Hereafter, we neglect the matter contribution because we are in deep radiation era. 
We can also neglect photon anisotropic stress since we focus on the epoch before 
recommbination. 

Eliminating g's from Eqs. ( |2.13| ) and ([2.14 ), we obtain following equation: 
3H 2 

-(i? 7 A 7 + R V A V + i? 7 A cx + R u tt u + i? 7 7r ex ) 



2$ 



k 2 



90/ 

+ — (Z-a) = 



(2.15) 



In this paper, we assume an adiabatic condition. In this case, 5P is expressed as 
5P = c 2 A, where c 2 is the adiabatic sound speed. Using the relations A 7 = 3<5P 7 , 
A u = 35 P v and A cx = 35P CX , Eq. (pT2| ) becomes: 



3V 2 

Z' + HZ + -— (i? 7 A 7 + R U A U + P 7 A C 
k 







From Eqs.( ^.l5| ) and ( |2.16| ), we obtain 



Z' + 3HZ + 2fc$ - 2Ua 



3U 2 
k 



(R^TT U + i? 7 7T, 



r 7 " ex ) 



Next we eliminate $ from Eqs. ( [2.11 ) and ( 2.17 ) using Eq.( [2.9|) : 

3H 2 



a' + 2Ua - krj 



k 



-{R v Tl v + -R 7 7T CX ) 



Z' + 3HZ -krj- (a' + 2Ho) 



3H 2 



(R U 7T U + -R 7 7T ex ) 



(2.16) 

(2.17) 

(2.18) 
(2.19) 
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From Eqs. (p.l8|) and (|2.19|) , we obtain 

Z' + 3HZ -2krj = 
Eqs. (|2.18|) and Q2.20|) together imply 



377" + 9Hrj' - kUo + k 2 V = -3U 2 {R v -k v + R^ir, 



•7 ' 1 ex / 5 



(2.20) 



(2.21) 



where we have used the relation 7/ = — 1(2 — a). After dividing Eq. (|2.21|) by R, we 



take the derivative of Eq.( [2.2lD with r: 



m- l r,'" - m- 2 H'i' + 9tj" + k 2 R~ x r{ - tfU^U'r) - ka' 

= -m'iRu-Ku + R^ cx ) - 3H(R v n' v + ii>ex) • (2-22) 



Finally we obtain a third order equation for rj from Eqs. (|2.18|) , ( |2.21| ) and (|2.22|) : 



3?ry" + (15 - 3n- 2 n')ri" + (isn + fc 2 ?rV + k 2 (i - n- 2 n')ri 

= {-m 2 - 3n')(R u ir u + # 7 7r ex ) - 3H(Ry v + # 7 < x ) . (2.23) 



In the radiation dominated epoch, R is approximately r Thus, we obtain the final 
result: 



k 2 r 2 2 
tV + 6rV + (6 + -T-W + -fcVr? 

= -2(R U 7T U + -R 7 7T CX ) - T(R U 7C' U + Ryir' ex ) 



(2.24) 



This is the equation for 77 sourced by the anisotropic stress. 

In order to solve the growth of the neutrino anisotropic stress, we need the 
relevant Boltzmann equation |1], ||, |20], |2l|j : 



< = k{\ qv -\Gf) + ^ka 



(2.25) 



where is the octupole moment. Eq. (|2.25| ) is valid after the neutrino decoupling. 
Eliminating a from Eq.( p.l8|) with Eq.( [2.25| ), we obtain the following equation: 



d 2 ir h 



+ 



dn,. 



+ -Ru^u + 0(kr) = --i? 7 7r e 



rf(lnr) 2 ' d\nr ' h^" v ' ~ v "'> 5 7 ex ' 
where the term of 0{kr) is negligible on super horizon scales. 



(2.26) 



3. Analytical Solutions 

To clarify the numerical evolution of the anisotropic stress we here describe an an- 
alytic model. We can set ir v ~ 7r£, ~ during the interval between the generation 
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Figure 1: (a) Numerical evolution of rj, R u ir u , i? 7 7r ex , and R u ir v + -R 7 7T ex . Neutrino 
decoupling occurs at t v ~ 9.3 x 10~ 5 Mpc. We set k = 10 Mpc , r T = 10 5 and 7r ex = 
1 ( r in < T ) and assume for convenience that the extra anisotropic stress was generated 
as a step function at r- ia = r u /r T . The neutrino anisotropic stress (red curve) cancels the 
extra anisotropic stress (violet curve). The curvature perturbation rj (green curve) stops 
growing after this cancellation. Panel (b) shows fit to the WMAP-5yr data of the primary 
scalar-mode and passive scalar mode power spectra induced by the extra anisotropic stress 
in the TT mode. This fit is compared with the standard ACDM model. The two power 
spectra are nearly indistinguishable. Note however that we adopted the spectral index 
which best fits the WMAP-5yr data for the extra anisotropic stress. 



epoch of the extra anisotropic stress, r in , and the epoch of neutrino decoupling, r v , 
i.e. T in <r<r u . 

First, we consider the case of r in < t < r v . Assuming n cx oc const., we find the 
following physical homogeneous and inhomogeneous solutions for Eq. ( |2.24| ), 

rj oc r _1 , const, (homogeneous solutions), 

n = — i? 7 7r ex lnr (inhomogeneous solution). (3.1) 

Although the homogeneous solution 77 oc r _1 is a decaying mode, we need this mode 
to calculate the growth of the curvature perturbations. Finally we obtain 

t](t) - r/(r in ) = -R y 7i ex (ln(r/r in ) + r in /r - 1) , (3.2) 

where we have used the initial condition rj = ^(r in ), 7/ = at r = r in . In this epoch, 
the curvature perturbation grows logarithmically. 

After neutrino decoupling, r > r v , the neutrino anisotropic stress must be 
evolved with the use of the Boltzmann equation. We can solve Eq. ([2.26 ) ana- 
lytically to obtain 

jr v = Cie (~l/2+V32^/5-Dlnr + ^(-1/2-^32^/5-1) Inr _ ? (3 3) 

where C\ and are constants of integration. Since the first two terms on the r.h.s. of 
Eq. (|3.3|) damp to zero, the neutrino anisotropic stress eventually cancels the extra 
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anisotropic stress asymptotically, 

Ry^v -> --R 7 7Tex • (3.4) 

Since r\ satisfies rj — — — a), we can rewrite Eq.( [2.25|) as follows: 



K = \v' + ^kZ + hq v -hGf (3.5) 
Assuming that kZ, kq v and kG v are negligible on super horizon scales, we have 

5 r 

v(t) ~ - g / <(r)df . (3.6) 
Thus, we obtain an asymptotic solution for 77 from Eqs. (|3.4|) and ( |3.6|) : 

17(1-) - r^) -|^r« • (3-7) 

After the complete cancellation of the anisotropic stress, therefore, the curvature 
perturbation stays constant. 

Finally, we obtain the result: 

T] -> T]pr + Tjpsi + 7] ps2 , 

r] psl = --R 7 7r ox (lnr r + r^ 1 - l) , 

?7 PS 2 = -5Ryir e J{8R„) , (3.8) 

where r T = r u /r in , r\ w is the primary perturbation, and ?7p S i represents the perturba- 
tions which are generated by 7r ex during the epoch of r in < r < r v . The quantity r] ps2 
is the perturbation induced by the growth of the anisotropic stress after neutrino 
decoupling. This provides insight into the numerical solutions described in the next 
section. 

4. Numerical Solutions 



We solved Eqs. fl2.24| ) and ( 2.26|) numerically. Fig. |l|(a) shows an example of the 



numerical evolution of 77 and the tt's. For purposes of illustration we in this figure 
have set k = 1(T 5 Mpc -1 , r T = 10 5 and 7r ex = 1 (r in < r). (Note that r T = 10 5 is 
not the same as the value adopted for figure 0(b), but is chosen here to make the 
graph easier to read). We have also assumed that the extra anisotropic stress was 
generated at r- m = r u /r T as a step function. This choice of a step function is a matter 
of convenience. In a realistic model, the appearance of this term would be a result 
of dynamics in the bulk dimensions. Nevertheless, as should be apparent from this 
plot, it does not particularly matter how the anisotropic stress grows in as long as 
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7r cx obtains its constant value sometime before neutrino decoupling. The numerical 
calculation agrees well with our analytical solutions [Eqs. ( p.4|) and (|3.8| )1. 

This asymptotic behavior to obtain constant rj is very similar to the standard adi- 
abatic mode whose r\ also remains constant on super horizon scales. These primary- 
like modes which are induced by the anisotropic stress, are called "passive modes", 
while the cancellation of the anisotropic stress generates "compensated modes" [|T^]. 

From the above discussions, we can conclude that if r T is large enough, the passive 
scalar mode can mimic the primary mode. If there is a primordial stress source, T^, 
which does not contribute to the background and satisfies = = 0, just like 
the dark radiation in the low-energy limit, the energy density and pressure grow as 
oc a -4 on superhorizon scales. Assuming that the anisotropic stress also scales as 
oc a -4 , we can calculate the CMB spectrum induced by an extra anisotropic stress 
as shown in Fig. |l](b). Here we set r T = 10 18 corresponding to the epoch at the end 
of inflation, and a/ | vr ex | 2 ~ 8.6 x 10~ 6 from a fit to the observed amplitude of initial 
curvature perturbation = 2.41 x 10~ 9 |I3. In addition, we adopt A ex ~ — 7r ex 



as in the case of a primordial magnetic field and we have used the best fit 
parameters from the WMAP-5yr analysis adopting the best fit spectral index for the 
extra anisotropic stress spectrum. 

As can be seen in Fig. |l|(b), the CMB anisotropies induced by the extra anisotropic 
stress is very similar to the primary spectrum induced by inflation. Note, however, 
that although we have generated a spectrum of the extra anisotropic stress to be 
consistent with the observation, this requires that the magnitude and spectrum be of 
the correct form to mimic the standard inflation-generated fluctuations. In general, 
the observed CMB spectrum could be a superposition of spectra generated by infla- 
tion (or some other mechanisms) and the anisotropic stress, as shown in Eq. (|3.8|) . 
Therefore, further study of the origin of the anisotropic stress and whether it can 
have the magnitude and spectrum consistent with the observed spectrum is highly 
desirable. 

The Gaussianity of the fluctuations in our model depends upon the source of 
the extra anisotropic stress. For example, a PMF spectrum is non-Gaussian. Since 
a normalized anisotropic stress a/ | 7r cx | 2 ~ 8.6 x 1CT 6 corresponds to about 9.3nG, 
which is the same order as the PMF estimated in Ref. [E3L one should take account of 
the passive components in scalar and tensor modes of the CMB. It might be expected, 
therefore, that the whole observed spectrum of the present CMB be explained by the 
passive scalar mode of the PMF. However, that is not the case because the tensor- 
to-scalar ratio of the passive PMF mode turns out to be greater than unity from 
previous work fl2"4| . Therefore a PMF is excluded by the WMAP-5 data. There is 
also another difficulty in that the spectrum of the PMF is not well known. 

We can understand the growth of curvature perturbations from the Einstein 
equation. The traceless part of the Einstein tensor is equal to the anisotropic stress 
in the universe, cf. Eq. (|2.2|). Since the traceless part of the Einstein tensor is written 
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in terms of the anisotropic expansion rate of the universe, <r, the anisotropic stress 
causes the universe to expand anisotropically. As the anisotropy of the universe 
increases, the curvature of the universe also increases. After neutrino decoupling, 
the neutrino anisotropic stress grows. This tends to decrease the anisotropic ex- 
pansion rate, Eq. ( |2.25| ), i.e. neutrinos cancel the extra anisotropic stress. After 
this compensation, the anisotropic expansion rate vanishes and the universe expands 
isotropically. In this phase, the curvature perturbation, r], is conserved. 



5. Conclusion 



We have found a simple solution for Eq. ( |2.24| ) and have shown that if there exists 
an extra anisotropic stress which scales as p 7 vr cx oc a -4 , the anisotropic stress from 
neutrinos tc u exactly cancels 7r ex . Before neutrino decoupling, the curvature pertur- 
bations grow logarithmically. After neutrino decoupling, they become constant on 
superhorizon scales just like the standard adiabatic mode of inflation. This is be- 
cause the total anisotropic stress vanishes via a cancellation. Thus, the resultant 
CMB spectrum is a superposition of the primary and passive modes. 

As an illustration we have considered the possibility that the passive scalar mode 
has a scale invariant spectrum. In this case the extra anisotropic stress might even 
produce a power spectrum similar to the observed CMB. This suggests a possible 
additional consideration in the determination of cosmological parameters from the 
CMB. Also we note that the Gaussianity of the CMB fluctuations depends upon the 
source for the extra anisotropic stress. Hence, should this extra anisotropic stress be 
present in the observed power spectrum, it might be detectable by non-Gaussianity 
in the fluctuations. Future CMB observations of non-Gaussianity may, therefore, 
help to constrain this possibility. 

To summarize, our purpose here has been to suggest that such a contribution 
to the observed spectrum may exist. Thus, further studies on the amplitude and 
spectrum of the extra anisotropic stress in brane-world cosmology are warranted. 
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